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Abstract. We review the expected science performance of the New 
Gravitational- Wave Observatory (NGO, a.k.a. eLISA), a mission under study 
by the European Space Agency for launch in the early 2020s. eLISA will sur- 
vey the low-frequency gravitational-wave sky (from 0.1 mHz to 1 Hz), detecting 
and characterizing a broad variety of systems and events throughout the Uni- 
verse, including the coalescences of massive black holes brought together by galaxy 
mergers; the inspirals of stellar-mass black holes and compact stars into central 
galactic black holes; several millions of ultracompact binaries, both detached and 
mass transferring, in the Galaxy; and possibly unforeseen sources such as the 
relic gravitational-wave radiation from the early Universe. eLISA's high signal- 
to-noise measurements will provide new insight into the structure and history of 
the Universe, and they will test general relativity in its strong-field dynamical 
regime. 

PACS numbers: 04.25.dg, 04.80.Nn, 95.55.Ym, 97.80.Af, 97.60.Lf, 98.35.Jk, 
98.62.Js, 98.80.Cq 

1. Introduction 

Over the last two decades, as many as 2,500 articles on space-based gravitational- 
wave (GW) detection included mentions of LISA (the Laser Interferometer Space 
Antenna) [HUMS], the space-based GW interferometer planned and developed together 
by NASA and ESA. This collaboration between the two agencies ended in early 2011 
for programmatic and budgetary reasons. In fact, LISA, as brought forth by the 
entirety of those papers, was more than a space project: it was the concept (and 
the cherished dream) of a space-based GW observatory that would explore the low- 
frequency GW sky, in a frequency band (10 _4 -1 Hz) populated by millions of sources 
in the Galaxy and beyond: compact Galactic binaries; coalescing massive black holes 
(MBHs) throughout the Universe; the captures of stellar remnants into MBHs; and 
possibly relic radiation from the early Universe. 

All along its evolution, the LISA design remained based on three architectural 
principles developed and refined since the 1970s: a triangular spacecraft formation 
with Mkm arms, in Earth-like orbit around the Sun; the continuous monitoring of 
inter-spacecraft distance oscillations by laser interferometry; drag-free control of the 
spacecraft around freely falling test masses, the reference endpoints for the distance 
measurements, achieved using micro-Newton thrusters. The current incarnation of 
this concept is eLISA (evolved LISA), a mission under consideration by ESA alone 
(under the official name of NGO, the New Gravitational- wave Observatory) for launch 
in 2022 within the Cosmic Vision program. 

The eLISA design would achieve a great part of the LISA science goals, as 
presented in [T], and endorsed by the 2010 U.S. astronomy and astrophysics decadal 
[4]. This article surveys eLISA's science performance (sensitivity, event rates, and 
parameter estimation), as scoped out by these authors in the spring and summer of 
2011, and as discussed in full in Ref. [5]- This article is organized as follows: in Sec. [2] 
we provide a very brief overview of eLISA and its GW sensitivity, while later sections 
are organized by science topics. In Sec. [31 we discuss the astrophysics of compact 
stellar-mass binaries in the Galaxy; in SecT^J the origin and evolution of the massive 
BHs found at the center of galaxies, as studied through their coalescence GWs; in Sec. 
[5j the dynamics and populations of galactic nuclei, as probed through the captures of 
stellar-mass objects into massive BHs; in Sec.|6j the fundamental theory of gravitation, 
including its behavior in the strong nonlinear regime, its possible deviations from 
general-relativistic predictions, and the nature of BHs; in Sec. [7j the (potentially new) 
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physics of the early Universe, and the measurement of cosmological parameters with 
GW events. Last, in Sec. [8] we draw our conclusions, and express a wish. 

2. The eLISA mission and sensitivity 

We refer the reader to [5] for a detailed description of the eLISA architecture. eLISA 
has a clear LISA heritage, with a few substantial differences. The eLISA arms 
will be shorter (1 Mkm), simplifying the tracking of distant spacecraft, alleviating 
requirements on lasers and optics, and reducing the mass of the propellant needed 
to reach the final spacecraft orbits. The orbits themselves may be slowly drifting 
away from Earth, again saving propellant, and the nominal mission duration will be 
two years, extendable to five. As much existing hardware as possible, including the 
spacecraft bus, will be incorporated from the LISA Pathfinder mission, scheduled for 
launch by ESA in 2014. The three spacecraft will consist of one "mother" and two 
simpler "daughters," with interferometric measurements along only two arms, for cost 
and weight savings that make launch possible with smaller rockets than LISA. (Note 
that LISA was to be built with laser links along the three arms, but it was not a 
requirement that they would operate throughout the mission.) 

The eLISA power-spectral-density requirement for the residual test-mass 
acceleration is S acc (f) = 2.13 x 10~ 29 (1 + 10~ 4 Hz//) m 2 s~ 4 Hz" 1 , while the position- 
noise requirement breaks up into S sn (f) — 5.25 x 1CP 23 m 2 Hz for shot noise, 
and 5 om „(/) = 6.28 x 1CT 23 m 2 Hz" 1 for all other measurement noises. With 
these requirements, eLISA achieves the equivalent-strain noise plotted in Fig. [IJ and 
approximated analytically by 

_ 20 4S acc (/) + S sn (/) + S onm (/) ( , f x\ 2 
5U) ~T 1? X [ 1+ \0A1c/2l)) ' (1 > 

where L = 1 Mkm, c is the speed of light, and S(f) has already been normalized to 
account for the sky-averaged eLISA response to GWs. This noise allows the detection 
of a strain of about 3.7 x 10~ 24 in a two-year measurement with an SNR of 1. The 
requirement on the useful measurement band is 10~ 4 Hz to 1 Hz, with a goal of 3 x 10~ 5 
Hz to 1 Hz. 



3. Compact binaries in the Galaxy 

The most numerous sources in the low-frequency GW sky observed by eLISA will 
be short-period binaries of two compact objects such as white dwarfs (WDs) or 
neutron stars (NSs). These systems have weak GW emission relative to the much 
heavier massive-BH binaries, but are numerous in the Galaxy and even in the Solar 
neighborhood. To date, astronomers have observed about 50 ultra-compact binaries 
with periods shorter than one hour, comprising both detached systems and interacting 
binaries where mass is being transferred from one star to the other. Wide-field and 
synoptic surveys such as SDSS and PTF (and in the future, PanSTARRS, EGAPS, and 
LSST) will continue to enlarge this sample [3 [8]. Interacting ultra-compact binaries 
with NS accretors are found by all-sky X-ray monitors and in dedicated surveys [3]. 

A large subset of known systems will be guaranteed verification sources for 
eLISA [TO]; their well-modeled GW signals will be detected within the first few 
weeks to months of operation, verifying instrument performance. The most promising 
verification binaries are the shortest-known-period interacting systems HM Cnc (with 
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Figure 1. cLISA equivalent-strain noise, averaged over source sky location and 
polarization, as a function of frequency. The solid red curve was obtained with 
the LISACode 2.0 simulator [B], while the dashed blue curve is plotted from Eq. 
pj. For comparison, the dotted green curve shows the LISA sensitivity. 



a period of 5.4 m [H]), V407 Vul (P = 9.5 m), and ES Cet and the recently discovered 
detached system SDSS J0651+28 (P = 12 m; [12]). 

eLISA will individually detect and determine the periods of several thousand 
currently unknown compact binaries (in our estimate, 3,500-4,100 systems 
for a two-year observation; [5] Q3]), while the combined signals of tens of millions 
unresolvable systems will form a stochastic GW foreground at frequencies below a 
few mHz ([15j[T6]; see Fig. [2]) About ~ 500 close or high-frequency (> 10 mHz) sources 
will be seen with large SNRs, allowing the determination of sky position to better than 
10 deg 2 , of frequency derivative to 10%, of inclination to 10 deg, and of distance to 
10%. This large sample will allow a detailed study of the Galactic population, which 
is poorly constrained by EM observations and theoretical predictions [17) . 

Detections will be dominated by double WD binaries with the shortest periods 
(5-10 minutes). Their mergers are candidate progenitors for many interesting systems: 
type la [T5] and peculiar supernovae [111120]; single subdwarf O and B stars, R Corona 
Borealis stars and maybe all massive WDs [5T] ; and possibly the rapidly spinning NSs 
observed as ms radio pulsars and magnetars [25]. These binaries are short lived, very 
faint for telescopes, and scarce (few thousand in the whole Galaxy), so GWs will 
provide a unique window on their physics. eLISA will determine their merger rate, 
constrain their formation, and illuminate the preceding phases of binary evolution, 
most notably the common-envelope phase. 

Common-envelope evolution is crucial to most binary systems that produce 
high-energy phenomena such as 7-ray bursts and X-ray emission, but our 
understanding of its physics and outcome is limited [23J [24] and challenged by 
observations [551 US] • The standard scenario is as follows. Most stars in the Universe 
arc in binaries, and roughly half of binaries are formed at close enough separations that 
the stars will interact as they evolve into giants or supcrgiants. Following runaway 
mass transfer, the companion of the giant can end up inside the outer layers (the 
envelope) of the giant; dynamical friction reduces the velocity of the companion, 
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Figure 2. Main figure: power spectral density of the stochastic GW foreground 
from Galactic binaries, before (blue) and after (red) the subtraction of individually 
resolvable systems, which are plotted as green and red/blue dots (for detached and 
mass-transferring systems) . A few known verification binaries are shown as white 
dots. The solid/dashed black curves trace instrument noise alone/with confusion 
noise. Spectra are shown for the TDI observable "X" (see, e.g., [13|): subtraction is 
simulated for a two-year observation and threshold SNR = 7; resolvable systems 
are placed a factor SNR 2 above the combined instrument and confusion noise. 
Inset: time series of the residual foreground, which carries information about the 
number and distribution of binaries in the Galaxy. 

shrinking the orbit and transferring angular momentum and energy into the envelope; 
the envelope eventually becomes unbound, leading to a very compact binary consisting 
of the core of the giant and the original companion [2 7) . 

eLISA will also test dynamical interactions in globular clusters, which produce 
an overabundance of ultra-compact X-ray binaries consisting of a NS accreting 
material from a WD companion. The eLISA angular resolution will be sufficient 
to distinguish WD binaries in clusters, verifying whether they are also plentiful. 

The eLISA measurements of individual short-period binaries will provide a wealth 
of information on the physics of tidal interactions and the stability of mass transfer. 
For detached systems with little or no interaction, the evolution of the GW signal is 
dominated by gravitational radiation: 

ho,M^f^D-\ /ocA^/ 3 / 1173 , f = YJ> ( 2 ) 

where h is the GW strain, / the GW frequency, Ai — (m 1 TO2) 3 ^ 5 /(rni + m^) 1 / 5 is 
the chirp mass with mi, to 2 the individual masses, and D is the distance. Thus, 
measuring h, /, and / (which will be possible in 25% of systems) provides and D; 
measuring also / (which may be possible for a few high-SNR systems) tests secular 
effects from tidal and mass-transfer interactions. Short-term variations are not likely 
to prevent detection |28j . and the precision of / and / determination increases with 
the duration of the mission. 
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Tidal interactions are possible when at least one binary component does not 
corotate with the orbital motion, or when the orbit is eccentric. Their strength is 
unknown [29 , and has important consequences on the tidal heating (and possibly 
optical observability) of WD binaries, as well as the stability of mass transfer. This 
process begins after gravitational radiation shrinks detached binaries to sufficiently 
close orbits (with P ~ a few minutes) that one of the stars fills its Roche lobe and 
its material can leak to the companion. Mass transfer can be self-limiting, stable, or 
unstable, depending on the resulting evolution of the orbit and of the donor radius. 
Unstable transfer leads to mergers; stable systems (the interacting WD binaries known 
as AM CVn systems, as well as ultra-compact X-ray binaries) will be observed - and 
counted - by eLISA in the early stages of mass transfer [30]. Efficient tidal coupling 
can return angular momentum from the accreted material to the orbit [29[ I31[ 132] , 
slowing the inspiral and increasing the fraction of WD binaries that survive the onset 
of mass transfer from 0.2% to 20% [33] , 

The unresolved foreground from Galactic binaries will provide a light 
additional noise component for the detection of loud broadband signals (see Fig. [2| , 
but it also contains precious astrophysical information. Its overall level measures 
the total number of binaries (mostly double WDs); its spectral shape characterizes 
their history and evolution; and its yearly modulation JM], together with the distance 
determinations from many individual systems, constrains the distribution of sources 
in the different Galactic components. Thus eLISA will probe dynamical effects in the 
Galactic center, which may increase the number of tight binaries [35] : it will measure 
the poorly known scale height of the disk; and it will sample the population of the 
halo [35] [TS] , which hosts two anomalous AM CVn systems and which may have a 
rather different compact-binary population than the rest of the Galaxy. Furthermore, 
the eLISA measurements of orbital inclinations for individual binaries, compared with 
the overall angular momentum of the Galaxy, will provide hints on the formation of 
binaries from interstellar clouds. 

eLISA will also constrain the formation rate and numbers of NS binaries and 
ultra-compact stellar-mass BH binaries, throughout the Galaxy and without EM 
selection effects. These numbers are highly uncertain, but as many as several tens of 
systems may be detectable by eLISA [331 ISZ] > complementing the ground-based GW 
observations of these same systems in other galaxies (and at much shorter periods) . 

More generally, the astrophysical populations and parameters probed by eLISA 
will be different from, and complementary to, what can be deduced from EM 
observations. For instance, eLISA will be sensitive to binaries at the Galactic center 
and throughout the Galaxy, while GAIA [35] will be limited to the Solar neighborhood; 
GWs encode distances and orbital inclinations, while EM emission is sensitive to 
surface processes. Dedicated observing programs and public data releases will allow 
simultaneous and follow-up EM observations of binaries identified by eLISA. 

4. Massive black-hole binaries 

(See [5] for a much deeper review.) According to the accretion paradigm [39 ] l40l 141"] . 
supermassive BHs of I0 6 -I0 9 M Q power quasars — active galactic nuclei so luminous 
that they often outshine their galaxy host, which are detected over the entire cosmic 
time accessible to our telescopes. Quiet supermassive BHs are ubiquitous in our low- 
redshift Universe, where they are observed to have masses closely correlated with key 
properties of their galactic host (see [42], and refs. therein) leading to the notion that 
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galaxies and their nuclear MBHs form and evolve in symbiosis (see, e.g., [43 ] [44 | 145] ) . 

In the currently favored cosmological paradigm, regions of higher-density cold 
dark matter in the early Universe form self-gravitating halos, which grow through 
mergers with other halos and accretion of surrounding matter; baryons and MBHs are 
thought to follow a similar bottom-up hierarchical clustering process [451 H7J 03 [45] 
[50] . MBHs may be born as small seeds (1O 2 -1O 3 M0) from the core collapse of the 
first generation of "Pop III" stars formed from gas clouds in light halos at z ~ 15-20 
[5TJ|S0]; or as large seeds (1O 3 -1O 5 M ) from the collapse of very massive quasi-stars 
formed in much heavier halos at z ~ 10-15 [521 153] ; or by runaway collisions in 
star clusters [S3]; or again by direct gas collapse in mergers [55J (See [55J [57J and 
refs. therein). The seeds then evolve over cosmic time through intermittent, copious 
accretion and through mergers with other MBHs after the merger of their galaxies. 

The cosmic X-ray background from active MBHs at z < 3 suggests that 
radiatively efficient accretion played a large part in building up MBH mass [551 [551 [50] . 
so information about the initial mass distribution is not readily accessible in the local 
Universe. By contrast, eLISA will measure the masses of the original seeds from 
their merger events. Furthermore, it is unknown [61] whether accretion proceeds 
coherently from a geometrically thin, corotating disk [55] (which can spin MBHs up 
to the J/M 2 = 0.93-0.99 limit imposed by basic physics [63, 64.) or chaotically from 
randomly oriented episodes [55] (which typically result in smaller spins). eLISA's 
accurate measurements of MBH spins will provide evidence for either mechanism [66] . 

After a galactic merger, the central MBHs spiral inward, together with their 
bulge or disc, under the action of dynamical friction, and pair as a pc-scale Keplerian 
binary [67, 681 [70} [TT] ; MBH binaries are then thought to harden into gravitational- 
radiation-dominated systems by ejecting nearby stars (assuming a sufficient supply) 
[72l [731 [74] or by gas torques and flows in gas-rich environments [75J [75J [77] ; the final 
binary coalescence is the most luminous event in the Universe (albeit in GWs). BH 
mergers have been explored only recently by numerical relativity [78] , showing how the 
mass and spin of the final BH remnant arise from those of the binary components, and 
predicting remarkable physical phenomena such as large remnant recoils for peculiar 
spin configurations 79J. The predicted coalescence rate in the cLISA frequency band 
ranges from a handful up to few hundred events per year, depending on theoretical 
assumptions ([80j [HI [82j [SJ El [851 (Ml [87]). 

eLISA will be sensitive to GW signals from all three phases of MBH coalescence 
(inspiral, merger, and ring-down [55]). To assess the eLISA science performance in 
this area, after experimenting with different waveform families, we modeled these 
signals with the "PhenomC" phenomenological waveforms [55] , which stitch together 
post-Newtonian (PN) inspiral waves |90j with frequency-domain fits to numerically 
modeled late-inspiral and ringdown waves. (PhenomC waveforms were derived for 
binaries with aligned spins; we adapt them to the non-aligned case by projecting the 
orbital angular momentum and individual spins onto the angular momentum of the 
distorted BH after merger.) 

The first metric of performance is the detection SNR, angle-averaged over sky 
position and source orientation, which is plotted in Fig. || as a function of total 
rest mass and cosmological redshift (left panel) and as a function of total rest mass 
and mass ratio for binaries at z = 4 (right panel). eLISA covers almost all the 
mass-redshift parameter space of MBH astrophysics: any equal-mass binary with 
Mtot = 10 4 -10 7 Mq (the crucial "middleweight" range inaccessible to EM observations 
beyond the local Universe) can be detected (with SNR > 10) out to the highest 
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Figure 3. Left: constant-level contours of sky- and polarization-averaged SNR 
for equal-mass non-spinning binaries as a function of total rest mass Mtot and 
cosmological redshift z. The SNR includes inspiral, merger and ringdown. Right: 
SNR contours as a function of Mtot and mass ratio q = mi/m.2- 
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Figure 4. Left: distribution of expected SNR for MBH mergers as a function 
of z, computed from the SE/LE/SC/LC metacatalog (see main text). Right: 
likelihood for the mixing fraction T . for an individual realization of mixed model 
J'SE + (1 - J")LE with T = 0.45 (see main text). 



redshifts, while equal-mass binaries with M tot > 10 5 M & are seen in detail as strong 
signals (SNR > 100) out to z = 5. Binaries with M tot > 10 5 M Q and mass ratios < 10 
are seen with SNR > 20 out to z = 4. 

To evaluate expected SNRs in the context of realistic MBH populations, we 
consider four fiducial scenarios (SE, LE, SC, LC) where MBHs originally form from 
Small (~ 100 Mq) or Large seeds (~ 10 5 M ), and where they subsequently grow by 
Extended or Chaotic accretion. (See [HI] for details; here we enhance that analysis 
by including random spin-orbit misalignments up to 20 deg in E models [HI])- For 
each scenario we generate multiple catalogs of merger events, and join them in equal 
proportions into a single metacatalog. Figure [1] shows the resulting distribution of 
SNR with z: eLISA will detect sources with SNR > 10 out to z < 10, a limit imposed 
by masses of the expected binary population as a function of z. 

For the same metacatalog, Fig. [5] shows the expected accuracy of parameter 
determination, estimated using a Fisher-matrix approach based on PN inspiral 
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Figure 5. Parameter-estimation accuracy (relative frequency of fractional or 
absolute errors over SE/LE/SC/LC metacatalog) for primary and secondary 
redshifted MBH masses and dimensionless spins (mi and ni2, si/mi and 0.2/^12 j 
respectively), luminosity distance Dl and sky position Af2. 

waveforms with spin-induced precession, augmented with PhenomC merger-ringdown 
waveforms. eLISA can determine the redshifted component masses (m rc d s hift = 
(l + z) m rost ) to 0.1-1 %; the primary-MBH spins to 0.01-0.1; and the secondary-MBH 
spins to 0.1 in a fraction of systems. (Compare with EM MBH-mass uncertainties ~ 
15-200%, except for the Milky Way MBH, and with very large MBH-spin uncertainties 
from Ken iron line fits |93j.) The errors in Dl have a wider spread, from a few percent 
to virtual non-determination, while sky position Q is typically determined to 10-1000 
deg 2 . Compared to previous published estimates for LISA, the accuracy in determining 
both Dl and f2 is reduced for eLISA by having interferometric measurements only 
along two arms (although three arms were always a goal, not a requirement, for LISA). 

The next order of analysis is to combine multiple MBH-coalcscence observations, 
resulting in a catalog of binary /remnant parameters, into a single inference about 
the mechanisms of MBH formation and evolution throughout cosmic history. 
This problem was analyzed extensively by Sesana and colleagues 94 in the context 
of LISA. We repeated their analysis for eLISA, by generating 1,000 catalogs of 
detected mergers (over two years) for each of the four SE/LE/SC/LC scenarios, 
and comparing the relative likelihood p(A vs. B) = p(A\C)/[p(A\C) +p(B\C)] for each 
pair of scenarios (A, B), for C = A or B. We considered only detections with SNR > 8, 
and used spinless, restricted PN waveforms. Table [T] shows our results for a relative 
likelihood threshold 0.95: for instance, the first row on the left shows that if SE is 
true, it could be discriminated from LE and LC in 99% of realizations, but from SC 
only in 48% of realizations; the last row on the left shows that LC could not be ruled 
out in 2% of realizations when SE or SC are true, but in 22% of realizations when LE 
is true. This degeneracy between accretion mechanisms is an artifact of the spin-less 
assumption; including information about the spin of the final merged MBH (which 
can be measured in 30% of detections) provides essentially perfect discrimination. 

Last, because no theoretical model will exactly capture the "true" formation and 
evolution history of MBHs, we investigated eLISA's ability of measuring the mixing 
fraction < J- < 1 in a mixture model FA + (1 — F)B that produces coalescence 
events with probability T from scenario A, and 1 — JF from B. For instance, for the 
case .FSE + (1 — J r )LE with T = 0.45, T can be measured with an uncertainty of 
0.1 (see right panel of Fig. Q. Although highly idealized, this example shows the 
potential of eLISA's observations to constrain MBH astrophysics along their entire 
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Table 1. Model discrimination with eLISA MBH-binary observations. The 
upper-right half of each table shows the fraction of realizations in which the row 
model would be chosen over the column model with a likelihood threshold > 0.95, 
when the row model is true. The lower-left half of each table shows the fraction 
of realizations in which the row model cannot be ruled out against the column 
model when the column model is true. In the left table we consider only the 
measured masses and redshift for observed events; in the right table we include 
also the observed distribution of remnant spins. 
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cosmic history, in mass and redshift ranges inaccessible to EM astronomy. 

In closing this section, we note that eLISA may also detect coalescences of BHs 
with masses of 10 2 -10 4 M (intermediate-mass BHs, or IMBHs). These events do not 
result from hierarchical galaxy mergers, but they occur locally under the extreme 
conditions of star clusters. IMBHs may form in young clusters by way of mass 
segregation followed by runaway mergers [95, 96, 97, 98, 99 ; IMBH binaries may form 
in situ jlOOj , or after the collision of two clusters |101l 1102) . Although the evidence 
for IMBHs is tentative |103l 1104) , eLISA may observe as many as a few coalescences 
per year |101j out to a few Gpc [53] ; it may also detect stellar-mass BHs plunging into 
IMBHs in the local Universe [T05] . 

5. Extreme-mass-ratio inspirals and the astrophysics of dense stellar 
systems 

There is of course one galactic nucleus, our own, that can be studied and imaged 
in great detail [1061 HHH QUI EU EHl UTT]. The central few parsecs of the Milky 
Way host a dense, luminous star cluster centered around the extremely compact radio 
source Sgr A* . The increase in stellar velocities toward SgrA* indicates the presence of 
a (4 ± 0.4) x 10 6 M central dark mass |111) . while the highly eccentric, low-periapsis 
orbit of young star S2 requires a central-mass density > 10 13 M Q pc -3 |112j : a density 
> 10 13 M Q pc -3 is also inferred from the compactness of the radio source [113] . These 
limits provide compelling evidence that the dark point-mass at SgrA* is an MBH 

nnnnnnis]. 

Unfortunately, the nearest external galaxy is 100 times farther from Earth than 
SgrA*, and the nearest quasar is 100,000 times farther, so imaging other galactic 
centers is prohibitive. It will however become possible with eLISA. This is because 
MBHs are surrounded by a variety of stellar populations, including compact stellar 
remnants (stellar BHs, NSs, and WDs) that can reach very relativistic orbits around 
the MBH without being tidally disrupted |116) . The compact stars may plunge directly 
into the event horizon of the MBH; or they may spiral in gradually while emitting GWs. 
These latter systems, known as extreme-mass ratio inspirals (EMRIs), will produce 
signals detectable by eLISA for MBH masses of 10 4 -10 7 M . Stellar-mass BHs should 
be concentrated in cusps near MBHs [117[ 1118) l98l 1119) 1120) and generate stronger 
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GWs thanks to their relatively larger mass, so they will provide most detections. 

EMRIs are produced when compact stars in the inner 0.01 pc of galactic nuclei 
are repeatedly scattered by other stars into highly eccentric orbits where gravitational 
radiation takes over their evolution [116] ; resonant relaxation caused by long-term 
torques between orbits increases the rate of orbit diffusion |121L I122j . although 
relativistic precession can hinder this mechanism [123J. EMRIs can also be made from 
the tidal disruption of binaries that pass close to the MBH [124] , possibly ejecting the 
hypervelocity stars observed in our Galaxy (see, e.g., |125j 1; and from massive-star 
formation and rapid evolution in the MBH's accretion disk [126 . Different mechanisms 
will lead to different EMRI eccentricities and inclinations, evident in the GW signal 

The detection of even a few EMRIs will provide a completely new probe of dense 
stellar systems, characterizing the mechanisms that shape stellar dynamics in the 
galactic nuclei, and recovering information about the MBH, the compact object, and 
the EMRI orbit with unprecedented precision |116j . Especially coveted prizes will 
be accurate masses for 10 5 -10 7 M MBHs in small, non-active galaxies, which will 
shed light on galaxy-MBH correlations at the low-mass end; MBH spins, which will 
illuminate the mechanism of MBH growth by mergers and accretion (see Sec. [4]); 
as well as stellar-BH masses, which will provide insight on stellar formation in the 
extreme conditions of dense galactic nuclei. The key to measurement precision is 
the fact that the compact object behaves as a test particle in the background MBH 
geometry over hundreds of thousands of relativistic orbits in a year; the resulting GW 
radiation encodes the details of both the geometry and the orbit |127|. 11281 1129} 1130] . 

To assess the eLISA science performance on EMRIs, we model their very 
complicated signals |131j using the Barack-Cutler (BC) phenomenological waveforms 
|132j . which are not sufficiently accurate for detection, but capture the character and 
complexity of EMRI waveforms. We complement this analysis with more realistic 
Teukolsky-based (TB) waveforms obtained by solving the perturbative equations for 
the BH geometry in the presence of the inspiraling body |133| ; these have been 
tabulated for circular-equatorial orbits and for some values of MBH spin |130U134j . 

To evaluate expected EMRI detection horizons and detection rates, we perform 
a Monte Carlo over 500,000 realizations of the source parameters, taking MBH rest 
mass in [10 4 ,5 x 10 6 ] Mq with a uniform logM. distribution; MBH spin uniformly 
in [0,0.95]; compact-body mass of 10 M Q , representative of a stellar-mass BH; orbit 
eccentricity before the final plunge uniformly in [0.05,0.4]; and all orbital angles and 
phases with the appropriate uniform distributions on the circle or sphere, with an equal 
number of prograde and retrograde orbits. We take the poorly known EMRI formation 
rate to scale with MBH mass as 400Gyr -1 (M./3 x 10 6 M Q )-°- 19 [EH H3H Q32] , and 
we distribute systems uniformly in comoving volume. Our assumptions are consistent 
with the MBH mass function derived from the observed galaxy luminosity function 
using the M,-a relation, and excluding Sc-Sd galaxies [1381 H391 1134j . 

We further assume an observation time of two years, consider EMRIs in the last 
five years of their orbit [B4], and require a detection SNR = 20 [HH EH [142]. The 
left panel of Fig. [6] shows the resulting maximum horizon redshift for BC waveforms, 
as a function of MBH rest mass — that is, it shows the z at which an optimally oriented 
source with the most favorable MBH and orbit parameters (as found in the Monte 
Carlo) achieves the detection SNR. Thus, EMRIs in the eLISA range will be detectable 
as far z — 0.7. By contrast, EM observations of 10 4 -10 6 M Q MBHs are possible in 
the local Universe out to z ~ 0.1. The right panel plots the distribution of SNRs as a 
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log, (M./M„) redshift z 

Figure 6. Left: maximum horizon redshift vs. MBH rest mass, BC waveforms 
(red curve); averaged horizon redshift vs. MBH rest mass, TB waveforms with 
a, /M m = and 0.9. Assumptions are given in the main text; the maximum is 
computed as the highest z with SNR > 20 in a given mass bin. Right: maximum 
EMRI SNR vs. redshift, BC waveforms. 

function of z, which shows that nearby EMRIs in the local Universe will yield SNRs 
of many tens. 

For comparison, the left panel of Fig. [6] shows also the horizons computed with 
sky- and orientation-averaged SNRs, using TB waveforms from circular-equatorial 
orbits with MBH spins a./M. = and 0.9. The difference between the BC and TB 
curves is consistent with the effects of sky-averaging: SNRs for optimally oriented 
systems are expected to be 2.5 times higher than averaged SNRs. The a./M, = 0.9 
systems are favored because high MBH spin allows for orbits closer to the event horizon 
and higher GW frequencies, which shifts the peak eLISA sensitivity to higher masses. 

The resulting number of expected eLISA detections over two years is ~ 50, 
as evaluated with the BC-waveform Monte Carlo, and ~ 30/35/55 (for a./M. = 
0/0.5/0.9), as evaluated with TB-waveform sky-averaged horizons. The higher TB 
event rate is explained by the inclusion of eccentric systems, which radiate more energy 
in the eLISA band, and it should be more reliable because of the broad sampling 
of source parameters. Remember however that EMRI rates arc highly uncertain 
[116, 135, 136, 123 . Even with as few as 10 events, the slope of the MBH mass function 
in the 10 4 -10 6 M range can be determined to 0.3, the current level of observational 
uncertainty |143j . 

Because EMRI waveforms are such complex and sensitive functions of the 
source parameters, these will be estimated accurately whenever an EMRI is detected 
|140[ 1141] I142j . In particular, we expect to measure the MBH mass and spin, as well 
as the compact-body mass and eccentricity to better than a part in 10 3 |132j . As an 
example, Fig.[7]shows the posterior distributions of the best-determined parameters for 
a z — 0.55 source detected by eLISA with SNR = 25, as computed with the Markov 
Chain Monte Carlo algorithm of |144j ; for this source, the luminosity distance Dj, 
would be determined to 1%, and the sky location to 0.2 deg 2 . Even with relatively low 
SNR, parameter-estimation accuracy is excellent. In general, we find that the eLISA 
and LISA parameter-estimation performance is very similar for EMRIs detected with 
the same SNR (but of course different distances) , so the reader can refer to treatments 
for LISA in the literature [T32l 11351 \TM . 
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Figure 7. Posterior probability plot for source parameters (MBH rest mass M, , 
MBH spin a., compact-body mass m, and orbit eccentricity at plunge e), in the 
SNR = 25 detection of a 10 + 10 6 M EMRI at z = 0.55, with a./M. = 0.7 and 

^plunge ~ 0.25. 



6. Precision measurements of strong gravity 

Einstein's theory of gravity, general relativity (GR), has been tested rigorously in 
the Solar system and in binary pulsars [1471 1148] : these tests, however, probe only 
the weak-field regime where the characteristic perturbative parameter e = v 2 /c 2 ~ 
GM/(Rc 2 ) is very small, ~ 10 -6 -10~ 8 (here v is the velocity of gravitating bodies, 
M their mass, and R their separation). By contrast, eLISA's GW observations 
of coalescing MBHs (Sec. [l} and of EMRIs (Sec. ^ will allow us to confront GR 
with precision measurements of its dynamical, strong-field regime, and to verify that 
astrophysical BHs are really the Kerr mathematical solutions predicted by GR. 

Before considering the GR tests possible with each of these sources, we note 
that, by the second half of this decade, second-generation ground-based detectors are 
expected to routinely observe the coalescences of stellar-mass BHs and (possibly) of 
asymmetric systems such as a NS inspiraling into a 100 A/ Q BH. However, they will do 
so with 10-100 times lower SNRs than eLISA (for the brightest sources), and for up 
to 1,000 times fewer GW cycles; thus, eLISA will test our understanding of gravity in 
the most extreme conditions with a precision that is two orders of magnitude better 
than that achievable from the ground. (Although most of the references cited in the 
rest of this section were developed for LISA, their broad conclusions are applicable to 
sources detected with comparable SNRs by eLISA.) 

All three phases of MBH coalescence offer opportunities for precision 
measurements. The year- long inspiral signals can be examined for evidence of a 
massive graviton, resulting in a frequency-dependent phase shift of the waveform, 
improving current Solar-system bounds |1491 1150) ; they can yield stringent constraints 
on other theories with deviations from GR parametrized by a set of global parameters, 
such as massless and massive Brans-Dicke theories [151LI152] . theories with an evolving 
gravitational constant 153 , Lorentz- violating modifications of GR (154] ; last, various 
authors have considered testing inspiral waves for hypothetical, generic modifications 
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of their amplitude and phasing [155] 1156] 1157] 1158] . 

The merger of comparable-mass MBH binaries produces an enormously powerful 
GW burst, which eLISA will measure with SNR as high as a few hundred, even at 
cosmological distances. The MBH masses and spin can be determined with high 
accuracy from the inspiral waveform; given these physical parameters, numerical 
relativity can predict the shape of the merger waveform, as well as the mass and spin 
of the final remnant MBH [159] . and these can be compared directly with observations, 
providing an ideal test of pure GR in a highly dynamical, strong-field regime. 

The frequencies and damping times of the quasinormal modes (QNMs) in the final 
ringdown |160] are completely determined by the mass and the spin of the remnant, 
and therefore can be used to measure them [1611 1162] , while their relative amplitudes 
hold information about the pre-merger binary |163j . again providing a check of 
consistency between GR predictions for the phases of coalescence. Furthermore, 
the measurement of at least two QNMs |162j will test the Kerr-ness of the MBH 
|164j against exotic proposals such as boson stars and gravastars [1651 11661 11671 1168] . 
Modifications of GR that lead to different emission would also be apparent [1691 1170] . 

EMRIs are expected to be very clean astrophysical systems, except perhaps 
in few systems with strong interactions with the accretion disk |1711 11721 1173] . or 
with perturbations due to a second nearby MBH or star |1741 1175] . Over day-long 
timescales, EMRI orbits are essentially geodesies of the background geometry; on 
longer timescales, the loss of energy and angular momentum to GWs causes a slow 
change of the geodesic parameters. In the last few years of their evolution, as observed 
by eLISA, EMRI orbits are highly relativistic (R < 10 R,) and display extreme forms 
of periastron and orbital plane precession. Indeed, EMRI GWs encode all the mass and 
current multipoles of the MBH [12711176] , which for a Kerr BH are uniquely determined 
by the mass and spin alone (another manifestation of the "no-hair" theorem). For 
EMRIs with SNR = 30, eLISA will measure mass and spin to a part in 10 3 -10 4 , and 
the mass quadrupole moment Mi to a part in 10 2 -10 4 , thus testing the no- hair 
theorem directly [129] , See [1771 1178] for reviews of different ways to test the nature 
of astrophysical BHs. 

Other tests of the Kerr-ness of the central massive object have been proposed: for 
a boson star, the EMRI signal would not shut off after the last stable orbit [179] : for a 
gravastar, QNMs could be excited resonantly [168] ; for certain non-Kerr axisymmetric 
geometries, orbits could become ergodic or experience resonances [1801 I181j ; for 
"bumpy" BHs, orbits would again carry distinctive signatures [1271 11821 11831 [T51] . 
Modifications in EMRI GWs would also arise if the true theory of gravity is 
in fact different from GR, as are dynamical Chern-Simons theory [1851 1186] . 
scalar-tensor theories (with observable effects in NS-BH systems where the NS carries 
scalar charge [151] 1187] ). Randall-Sundrum-inspired braneworld models [188] 1189] . 
theories with axions that give rise to "floating orbits" [1901 1191] , as well as generic, 
phcnomenologically parametrized theories [192] . 

7. Cosmology and new physics from the early Universe 

GWs produced after the Big Bang form a fossil radiation: expansion prevents them 
from reaching thermal equilibrium with the other components because of the weakness 
of the gravitational interaction. Thus, relic GWs carry information about the first 
instants of the Universe. If their wavelength is set by the apparent horizon size 
c/H* = c(a/a)„ at the time of production, when the temperature of the Universe 
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Figure 8. (From I193| .) Spectra of stochastic backgrounds from cosmic strings 
for large loops (0.1 the horizon size, solid lines), for two values of the string tension 
G/x/c 4 spanning a range of scenarios motivated by braneworld inflation; and for 
small loops (with size a = 50eGfi, dashed line). The cosmic-string spectrum is 
distinguishably different from that of first-order phase transitions or any other 
predicted source: it has nearly constant energy per logarithmic frequency interval 
over many decades at high frequencies, and falls off after a peak at low frequencies, 
since large string loops are rare and radiate slowly. Cosmic strings may also 
produce distinctive bursts, produced by a sharply bent bits of string moving at 
nearly the speed of light \jl2[ [2131 [2141 |2T5| . 



is the redshifted frequency is 



/.^XIS.^P, (3) 

so the eLISA frequency band corresponds to the horizon at and beyond the Terascale 
frontier of fundamental physics. This allows eLISA to probe bulk motions at times 
about 3 x 10~ 18 -3 x 10 -10 s after the Big Bang, a period not directly accessible 
with any other technique. Taking a typical broad spectrum into account, eLISA has 
the sensitivity to detect cosmological backgrounds caused by new physics at energies 
~ 0.1-1000 TeV, if more than a (modest) fraction ~ 10~ 5 of the energy density is 
converted to GWs at the time of production. 

Various sources of cosmological GW backgrounds are presented in detail in 



[193 . They include first-order phase transitions, resulting in bubble nucleation and 
growth, and subsequent bubble collisions and turbulence |194[ I195[ I196( I197[ 1198] : 
the dynamics of stabilization for the extra dimensions required by superstring theory 
[19911200"] . which may also appear as non-Newtonian gravity in laboratory experiments 
at the sub-mm scale; networks of cosmic (super-)strings [2011 1202] . which continuously 
produce loops that decay into GWs (see Fig. [8]); the transition between inflation 
and the hot Big Bang in the process of preheating [2031 H2H HH US EOT]; and 
the amplification of quantum vacuum fluctuations in some unconventional versions of 
inflation [2081 12091 [2TD] . Although the two-arm eLISA does not provide a Sagnac TDI 
combination [211) to calibrate instrument noise against possible GW backgrounds, 
the clear spectral dependence predicted for some of these phenomena provides an 
observational handle, as long as the background lies above the eLISA sensitivity curve. 



Low-frequency gravitational-wave science with eLISA/NGO 



1G 



As discussed in Sec. |4j observations of GWs from MBH binaries probe the 
assembly of cosmic structures. In addition, binaries can serve as standard sirens 
to measure cosmological parameters [2161 1217] because, as discussed around Eq. 
([2]), measuring the amplitude and frequency evolution of a binary signal yields the 
absolute luminosity distance to the source. However, binary GWs cannot provide 
the source's redshift unless the other source parameters are known independently 
(because the rest mass of the binary is the only length/time scale in the waveform, 
the frequency evolution of a redshifted signal is indistinguishable from the signal from 
a heavier binary). The optical redshift of the host galaxy can be obtained if an EM 
counterpart to MBH coalescence is observed (see, e.g., [2181 12191 1220] . and |221j for a 
recent review). 

While there are many uncertainties in the nature and strength of such 
counterparts, some may be observable in the local Universe. At z < 1, we expect that 
eLISA MBH-inspiral measurements could provide sky locations to better than 400 deg 2 
for 50% of sources, and to 10 deg 2 for 11%. (The inclusion of merger and ringdown in 
the analysis should further improve these numbers.) Such large areas will be covered 
frequently and deeply by optical and radio surveys such as LSST |222j and the VAST 
project |223j . identifying sufficiently distinctive transients. The accurate knowledge 
of the counterpart's redshift and position would then improve the uncertainty of 
GW-determined parameters, with Dl known to 1% for 60% of sources, and 5% for 
87%. Such precise luminosity distance-redshift measurements will be complementary 
to other cosmographical campaigns |224[ [225], and will improve the estimation of 
cosmological parameters. Even without counterparts, one may proceed by considering 
all possible hosts in a distance-position error box, and enforcing consistency between 
multiple GW events [US]; this should be possible for MBH binaries (and EMRIs [227] ) 
in the local Universe, yielding the Hubble constant to a few percent. 

8. Conclusions 

While LISA was always meant to be the definitive mission in its frequency band, eLISA 
is being designed to provide the maximum science within a cost cap. Nevertheless, 
as described above, eLISA will achieve a great part of the LISA science goals. It will 
represent the culmination of twenty years of exciting, painstaking work, pioneering 
the new science of observational low-frequency GW astronomy. It will truly begin to 
unveil the hidden, distant Universe. May it fly soon, and safe. 
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